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INTRODUCTION 

This report is on heat-pipe research carried out for Ames Research 
Center under contract number NAS 2-8310. The wo~k focuses on the mathe- 
matical modeling of three critical mechanisms of heat-pipe operation. 

In particular. Section 1 is on the effect that excess liquid has on heat- 
pipe performance; Section 2 is concerned with the calculation of the 
dryout limit of circumferential grooves; and Section 3 is on an efficient 
mathematical model for the calculation of the viscous -inertial interaction 
in the vapor flow. These mathematical models are incorporated in the com- 
puter program GRADE II, which is described in Reference (4). 
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1.0 EXCESS LIQUID IN HEAT-PIPE VAPOR SPACES 

There are several reasons why heat pipes often operate with more liquid 
than necessary to just saturate the wick structure. For example, when a 
heat pipe is designed to operate over a wide temperature range, it is neces- 
sary to calculate the fluid charge for the temperature extreme that requires 
the most fluid. This usually occurs at the lowest operating temperature due 
to the increase of liquid density with a decrease in temperature; however, 
it may occur at the highest operating temperature due to the increase of 
vapor density with an increase in temperature. In either case, for operating 
temperatures other than the extreme for which the charge just saturates the 
wick, excess liquid will be present. This source of excess liquid can be 
particularly significant for gas-loaded variable-conductance heat pipes for 
use with spacecraft radiators. A relatively high fluid charge is required 
when the gas blocks a large fraction of the total heat-pipe length and the 
liquid in the gas-blocked region falls very low in temperature. 

Even when the entire heat pipe operates in a narrow temperature range, 
excess liquid is often intentionally introduced as a safety factor. For 
some wick types such as axial grooves, a small undercharge can significantly 
degrade the transport capacity while a small overcharge has little adverse 
effect. It is particularly important to ensure sufficient liquid in arterial 
or tunnel -wick heat pipes because even a slight undercharge makes priming 
impossible. 

Excess liquid affects heat-pipe performance in various ways. It resides 
in the heat-pipe vapor spaces and forms liquid fillets along corners and a 
puddle along the bottom. These fillets and puddles provide low-permeability 
flow paths in parallel with the wick structure. Their size at a given axial 
location depends on the local value of the vapor-liquid pressure difference, 
or capillary stress. Thus, excess liquid usually has the greatest effect in 
the condenser region where the stress is low. 

In addition to contributing to the axial transport of liquid, a puddle 
reduces the hydrostatic stress that the wick must support. Although excess 
liquid increases axial -transport capacity, it degrades the heat transfer, 
because fillets and puddles block condensation surface. The condenser end 
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of a vapor space can even completely fill a liquid slug that blocks con- 
densation heat transfer in that region altogether. An accurate mathematical 
model of excess-liquid behavior is particularly needed for the graded- 
porosity-wick heat pipe, where to maximize its capacity the wick is begun 
at the condenser end with the highest possible porosity. If the porosity 
is too high, however, the wick will not generate a sufficient capillary 
pressure to prevent a liquid slug from forming. 

The mathematical modeling of excess liquid is complicated by the fact 
that for typical heat pipes one cannot neglect variations of hydrostatic 
pressure with position on the free surface of the liquid. An earlier 
mathematical model in the heat-pipe computer program MULTIWICK^^^ neglected 
these variations. Although the resulting simplification allowed application 
quite general heat-pipe cross-sectional geometries, at low values of capil- 
lary stress the results are only qualitatively correct. To retain the 
hydrostatic pressure variations in the present work, it is necessary to 
select specific types of cross-sectional geometries for the vapor spaces. 

We have selected the circular vapor space which occurs in axially-grooved 
heat pipes and simple heat pipes with wick-lined wall, and "Dee-shaped" 
vapor spaces formed by slab wicks in circumferentially grooved tubes. 

Section 1.1 presents methods for calculating the contribution of excess 
liquid to axial transport, and Section 2.2 presents both theoretical and 
experimental results for conditions under which liquid slugs form in the 
vapor spaces. These results also apply to arteries, so we digress in Section 
1.2 to discuss criterion for the priming of open-and-closed arteries. 

1.1 FILLETS AND PUDDLES IN VAPOR SPACES 

The three geometries considered are shown in Figure 1-1, where typical 
fillets and puddles are drawn as they might appear at a specific axial 
location of a heat pipe operating in earth gravity. The key parameter that 
governs the configuration of the excess liquid is vapor-liquid pressure 
difference. Since it varies hydrostatically transverse to the heat pipe, 
in order to have a unique value at every axial location, we take its value 
along the centerline of the heat pipe anc define it as the capillary stress. 
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As the stress decreases, the size of the fillets and puddles increases. If 
the stress is sufficiently low, the vapor space will bridge and a liquid 
slug will form, which is the subject of the next section. Vie now present 
methods for the calculation of the free-surface shape of a fillet or puddle, 
its cross-sectional area, and its free and wette * perimeter. 

The equation governing the free surface of a liquid under the action 
of surface tension is 

aP = a(l/R^ + I/R 2 ), (1.1) 

where aP is the local pressure difference between the gas and liquid, ois 
the surface tension, and R^ and R 2 are the local radii of curvature in any 
two orthogonal planes that contain the normal to the surface. As the first 
of these, we take the x-y plane of the heat-pipe cross section. By assuming 
that the cross-sectional shape of the fillet or puddle changes gradually 
with distance along the heat pipe, we can neglect the reciprocal of the 
second radius of curvature compared to that of the first, and Eq, (1.1) 
reduces to 


aP = o/R 


( 1 . 2 ) 


A differential equation is obtained for the free-surface shape by 
substituting for aP the hydrostatic pressure variation with y and by sub- 
stituting for R the Cartesian formula for the radius of curvature, which 
results in 


pgy 




[1 + (dy/dx)^]^/^ 


(1.3) 


Here p is the difference in the liquid and vapor densities, and g is the 
gravitational acceleration. Also, the location of the coordinates is 
chosen such that the gas-liquid pressure difference is zero when y = 0. 
Equation (1.3), however, presents difficulties because one cannot integrate 
readily through regions where dy/dx is infinite. To remedy this, we use 
the transformation 


dy/dx = tan a (1 .4) 
to transform Eq. (1.3) into two differential equations where the dependent 
variables are x and y and the independent variable is the inclination a of 
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the free-surface curve. The resulting equations are 



2 

2y dx/do = a cos a. 

(1.5) 


2 

2y dy/dot = a sin a. 

(1.6) 

where 

a = /2a 7p 

(1.7) 


is the capillary constant. Equation (1.6) is readily integrated, which 
results in 

y = /yQ - a (cos a - cos Oq) (1.8) 

where the curve has an inclination at a point (x^, y^). Witn the 
solution (1.8), we can eliminate y from Eq. (1.5), which gives 

dx/da = a^ cos a (1.9) 

2 - a^(cos a- cos oIq) 

This equation can be integrated in terms of elliptic integrals by first 
using the transformation = 2'}>-if, which allows us to rewrite Eq. (1.9) as 

2 

dx = ^a ] d (1.10) 

/I - r sinS 


where 


k = 


2a^ 

p 2 

a (1 + cos a„) + y„ 

0 0 


( 1 . 11 ) 


One can easily verify that the solution to Eq. (1.10) in terms of elliptic 
integrals of the first kind F(k,(|)) and of the second kind E(k,(j>) is 


X = ( v^a/k) {(l-k2/2)[F(<^.k)-F(^^,k)] 

- [E((|>,k) - E((j.Q.k)]} + Xq (1.12) 

where is the value of <{> when a = a^. 

Equations (1.8) and (1.12) define a single free-surface curve that 
passes through the point (x^^, y^) at an inclination a^. The task now is to 
fit a segment of that curve inside the cross-sectional geometry of the heat 
pipe such that it represents a fillet or puddle. As seen in Figure 1-1, in 
the case of a fillet the boundary conditions require that, for a zero wetting 
angle, one end of the segment is tangent to the slab wick and the other is 
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tangent to the tube wall. We will describe how the solution is found for 
the puddle and for one particular fillet. The method for the other fillets 
is the same except for the details. 

Consider first, as shown in Figure 1-2, the solution for the puddle. 

The point at the bottom of the puddle is taken as the point (x^, y^), where 
Xq = 0, = 0, and y^, which sets the gas-liquid pressure difference at 

the surface of the puddle, is a selected positive value. Points (x,y) 
along the free-surface curve which are calculated from Eqs. (1.8) and (1.12) 
for increasing values of a, are checked to see if the condition of tangency 
with the tube wall is met, which is 

X - R.J. sin a = 0 (1.13) 

Figure l-2(b) shows that there is a second tangency condition for which 
Eq. (1.13) is satisfied. Thus the same free-surface curve can represent 
two puddles. 


y y 



Figure 1-2. Tangency Condition: x = Rj sin a 
(a) a < u/2 (b) a > •tr/2 
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The capillary stress (gas-liquid pressure difference at the center of 
the tube) is given by 

pgy^ = pg(y + Rj cosa) (i.i4) 

where y and a are the values at the point of tangency. A series of puddles 
are calculated for different values of capillary stress by starting with 
different values for y^. 

The bottom fillet that forms in a heat pipe with a vertical slab wick 
[Figure 1-1 (c)] is calculated exactly the same as a puddle with Eqs. (1.13) 
and (1.14) applying except that the initial condition is x^ = T/2 and 

= - a/2. Again two fillets are possible from a given free-surface curve. 
The other types of fillets depicted in Figure 1-1 are calculated by a simi- 
lar method except for these only a single fillet is possible for a given 
free-surface curve. 

We now turn to the calculation of areas, wetted perimeter and free 
perimeter. Although the free perimeter is not used in the calculation of 
the effective permeability of the excess liquid, it is an important para- 
meter for future work on vapor-liquid interaction. An elemental length 
dS^ of the free perimeter is given by 

dS^ = v^x^ + dy^ (1.15) 

With the use of Eqs. (1.5), (1.6) and (1.8), we rewrite Eq. (1.15) as 

a^dg 

^^f " 2 ^ 0 ^ - a^(cosa - cos g^) (1.16) 

As in the solution of Eq. (1.9), we use the transformation a = 2<t>-Tr to 
write dS^ as 

dS, = »^ka / “ — (1.17) 

/i - k*^ sin <t) 

where k is given by Eq. (1.10). Equation (1.16) integrates to 

= .^ ka [F( 4 ,,k) - F(<|.^,k)] (1.18) 
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For the calculation cf the area, refer to Figure 1-3 (a). Although a 
particular fillet is shown, the expressions derived below apply directly 
to the puddle and other fillets as well. The area desired is 
where the subscripts denote points and lines bounding the area. It is 
more convenient, however, to calculate A^,o* Then, the desired area is 

\ " ^023 * ^^2b3 ■ ^23^ ' ^oa2 (1.19) 



Figure 1-3. Calculation of the fillet area A^. 

(a) Diagram defining points, lines and vectors. 

(b) A blow-^p of Figure 1-3 (a) showing the 
vector q that connects (x , y ) to a point 
(x,y) on the free-surface curve. 


The area of triangles Aq 23 and A ^23 
coordinates for their vertices: 


are calculated in terms of the 


^23 = 7 I^V3 - 


( 1 . 20 ) 


^023 = j |(x2y3 - y2X3) - (Xpy3 - yQX3) + - yoX2)| (1.21) 
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The area of the sector is 1/2 e Rj*^, where o is the anqie between 

vectors ®nd : 

► 'f 

0 = cos"^ '’31 • *^21 


IP31I IP21 


^ (X3 - x^)(x2 - x^) + (y 3 - y^)(y2 - y^ ) 

/[(X 3 - x^)^ + (y^ - yil^][(x 2 - x^)^ + (y^ - y^)^” (1.22) 

It remains to calculate Cotuider Figure 1-3 (b) where the vector q 

connects the point (xQ.y^) to a general point along the free-surface curve. 
The area between q and the curve is A. We now calculate the increase in 
area dA when q is incremented an amount corresponding to the increment du. 
The incremental area dA is the area of the triang’c bounded by q and dq, 
thus 


dA = ^1 q X dq 


With q = (x-x^) T + (y-y^) J and dq = dx t + dy j, we find 


dA = ^|dy(x - Xq) - dx(y - y^)| 


(1.23) 


(1.24) 


Equations (1.5) and (1.6) are used to obtain 

^ " ly ■ ’^o^ ^ - i'J - yo) cos ot] (1.25) 

Since x is given in terms of elliptic integrals, there is no hope of 
obtaining an analytic integral for A. Therefore, the integration is 
carried out numerically from to the point where the boundary condition 
on the tube wall is met. 

Referring to Figure 1-3 (a), we see that the wetted perimeter is given 
simply by = y^ - y 3 + R.J. e. (1.26) 

(0) 

The computer FILLET' ' was written to carry out the calculations 
outlined in this section. Since we could not find a general-purpose sub- 
routine for elliptic integrals, Eqs. (1.9) for x and (1.17) for S^ are 
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integrated numerically along with Eq. (1.24) for A. The output from 
FILLET is a table of effective hydraulic diameter and the total area 
of fillets and a puddle for the range of capillary stress that can exist 
in a heat pipe. 

To illustrate the calculations, we focus attention on the puddle in 
a circular vapor space. Two parameters are required to fix the puddle 
configuration. The first is the Bond number 

B = pgR^/o, (1.27) 

which is the ratio of hydrostatic forces to surface tension. The second 
paranteter is the dimensionless capillary stress 

S=Po/p9Rt, (1.28) 

where is the gas-liquid pressure difference at the center of the tube. 
For a geometrical interpretation of the dimensionless stress, imagine that 
it is due solely to hydrostatic pressure = pgh. Thus, the stress is 
simply the ratio of the height h of the tube center above the flat surface 
of a liquid reservoir to the tube radius Rj. When the stress is negative, 
the tube center lies below the reservoir surface. For the case of fillets 
in a heat pipe with a slab wick, a third parameter, the ratio of the wick 
thickness T to the tube radius Rj, is required. 

The results for a puddle are displayed in Figure 1.4 (a) - (d) for 
values of the Bond number of 0.1, 1., 10., and 100. For each Bond number, 
the puddle configuration, its area, free perimeter and wetted perimeter 
are shown for a range of stress. The puddle area is made dimensionless 
by dividing by the c»'oss-sectional area of the tube, and the free and 
wetted perimeter are made dimensionless by dividing by the circumference. 
Thus we see that the dimensionless area and wetted perimeter approaches 
unity as the puddle nearly fills the tube. For high values of stress when 
the puddle is shallow, the free perimeter and wetted perimeter are nearly 
equal . 
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(e) Bo = 10 



0.00 0-*0 0.80 1-20 
DIMENSIONLESS STRESS VI /R 




.Jt \ 


(d) Bo= 100. 



_0 ,40 o'.OO 0.40 0-80 1 -20 

DIMENSIONLESS STRESS Vl/R 





Figure 1-4. Fillet configurations in horizontal circular cylinders. 

Dimensionless fillet area X, wetted perimeters^, and fr 
perimeter Tr are plotted as a function of dimensionless 
stress V.|/R. An asterisk denotes the neutrally stable 
meniscus! 


POLDOUT JltAME 
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As will be discussed in the next section, those liquid configurations 
for which the fillet area increases with increasing stress are unstable 
and cannot exist. Note that this is the case for all puddle configurations 
when B = 0.1. The critical-free surfaces, which are neutrally stable are 
denoted in the figure. For stresses below the critical stress, a liquid 
slug forms. 

1.2 EXCESS LIQUID SLUGGING OF VAPOR SPACES AND THE PRIMING OF ARTERIES 

In this section we focus on the conditions under which a liquid slug 
forms in a vapor space, which, as mentioned previously, also applies to 
the priming of arteries. There are two distinct situations to consider; 
vapor spaces with and without solid end boundaries. The former occurs in 
most simple heat pipes, such as an axially grooved or a slab-wick heat 
pipe that have end caps. The latter occurs primarily in gas-loaded heat 
pipes where the condenser end protrudes into a gas reservoir and is other- 
wise uncapped. We consider this situation first. 

1.2.1 Slugging of Open Vapor Spaces 

Consider an experiment, as shown in Figure 1.5, where a horizontal 
open-ended porous tube is lowered into a liquid pool. The capillary stress 
is measured by the elevation h of the tube center above the pool. If the 
Bond Number is sufficiently large, a puddle of cross-sectional area A will 
form before the tube slugs with liquid. As h is decreased A increases. 

The tube will abruptly slug with liquid when dA/dh = <=°. This slugging 
criterion, shown in Figure 1-6 by the dashed line, is found from the cal- 
culations of the program FILLET such as those displayed in Figure 1-4. 

In Figure 1-6, instead of plotting the critical height to the center of 
the tube, we plotted the critical height to the top, which eliminates 
negative values and allows the use of logarithmic graph paper. 

If the Bond number is sufficiently small, then dA/dh is positive for 
all puddles, and hence none are stable. Therefore, slugging abruptly occurs 
as soon as a puddle begins to form on the bottom of the tube. Such a puddle 
has a radius of curvature nearly equal to the tube radius and thus, by 
Eq. (1.1) the gas-liquid pressure difference at the tube bottom is a/Rj: 
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pg (h - Rj) = a/Rj. 

or in dimensionless form 

h/Ry - 1 = 1/B. (1.29) 



Figure 1-5. When the cylindrical tube is lowered sufficiently, a liquid 
slug will abruptly form. 


- 13 - 



26263-6025-RU-00 



Figure 1-6. Slugging criterion for open-ended cylindrical tubes. 
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This criterion for slugging is displayed in Figure 1-6 by the solid line. 
The intersection of the dashed and solid lines at a Bond number of .33 
represents the point at which a stable puddle is just possible. The ana- 
lytical calculation of this critical Bond number poses an intriguing 
academic challenge. Also displayed in Figure 1-6 is data taken by Johnson 
and Saaski^^^ on the priming of open-ended glass and screen walled tubes 
ranging from 0.51 mm to 3.2 mm with 2-propanol and Freon 113 as the liquid. 

The condition for open-ended slugging of "Dee-shaped" vapor spaces 
can also be found by the program FILLET, which calculates the fillet areas 
as a function of stress S. Slugging occurs when one of the two fillets, 
or the puddle if one is possible, becomes unstable, that is dA^/dS = ® . 


1.2.2 Slugging in a Closed Vapor Space 

In a vapor space with a closed end, a slug first begins to form at 
the end, and then grows progressively longer. Such a slug is depicted 
in Figure 1-7. The closed-vapor-space slugging criterion is more import- 
ant than the open criterion, therefore, in this section we present detailed 
results for the "Dee" shape as well as the circular vapor space. 

The critical capillary stress for the formation of a slug is higher for 
the closed vapor space owing to the fact that the free surface develops 
curvature in two orthogonal planes rather than just one. In fact, if the 
stress is sufficiently low for a liquid slug to form in an open vapor space, 
then to clear the slug the stress must be increased to the critical stress 
for slugging of the closed vapor space. For example, in the design of a 
graded-porosity wick, even if the condenser end opens into a gas reservoir, 
the slugging criterion for a closed vapor space should be used, because the 
condenser-end porosity must be low enough to generate a sufficient capillary 
pressure to clear an existing slug. Similarly, while the maximum stress to 
initiate priming of an open artery is set by the slugging criterion for open 
vapor spaces, the capillary pumping, that is, the stress necessary to initiate 
emptying, is set by the slugging criterion for closed vapor spaces. 
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The mathematical modeling of the slug that forms in the closed vapor 
space is difficult because the free surface varies in three dimensions. 

We consider first the circular vaoor space and employ an approximation 
suggested by Johnson and Saaski^^' to reduce the problem to a two- 
dimensional one. Applying Eq. (1.1) to the free surface shown in Figure 
1.7, we take as the local radius of curvature in the y-z plane which 
bisects the tube vertical Iv and R 2 as the local radius of curvature in 
plane orthogonal to the y-z plane that contains the local normal to the 
surface. The approximation is to use a constant characteristic value for 
R^, aiid for lack of a better candidate this constant value is taken as 
the tube radius Rj. The equation governing the meniscus shape in the z-y 
plane is developed exactly the same as Eqs. (1.5) and (1.6), except the 
constant term 1/Rj is carried along. Thus, instead of Tq. (1.6), we obtain 

(2y - l/Rj) ^ = a2 sin a (1.30) 

Referring to Figure 1.7, we see that the boundary conditions for the free 
surface curve is 

a = oaty = h-Rj + d 

and (1.31) 

a = -TTaty = h + Ry 

In case slugging occurs without a puddle forming first, we take d = 0 . 
Integration of Eq, (1.30) between the above boundary conditions results in 

d2 + d[2h - 2Rj - a2/Rj] + 4[a2 - R^j-h] = 0 (1.32) 

The no-puddle slugging condition is (d = 0 ) 

h/R-,- = 2/B (1.33) 

Equation (1.29) gives the condition for a puddle just about to form. The 
two curves given by Eqs. (1.29) and (1.33) intersect at B = 1. Thus, for 
B < 1, the slug forms without a puddle and Eq. (1.33) is the appropriate 
criterion. For B > 1 , we use a modification of the program FILLET to 
calculate the puddle depth as a function of h/Rj and B: 
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d/Rj = f(h/R^. B) (1.34) 

This relation is used to eliminate d Eq. (1.32). The resulting slug- 
ging criterion for B > 1 alo.ig with the criterion for B < 1 given by Eq. 
(1.33) are displayed in Figure 1.8. 



Figure 1.7. Liquid Slug with a Puddle of Depth d on 
the Bottom of the Tube 

The validity of the approximation of Johnson and Saaski can only be 
assessed by comparing the resulting slugging criterion to experimental 
results or to an exact numerical solution. We carried out an experiment 
with a series of glass tubes that varied in diameter from 2 to 30 mm. As 
shown in the Figure 1.9, one end of each tube was drawn down to half the 
original diameter and the tapered section was bent 90°. The set of tubes 
were placed on a level platform inside a large glass jar with the tapered 
sections down. The jar was slowly filled with acetone. At a critical 
level, a liquid slug, which initiated in the tapered section, moves ra- 
pidly through the smallest tube. The height of the critical level is 
measured with a cathetometer. The jar is slowly filled again until the 
slug moves through the second tube, and the height is measured again, and 
so forth. From these data, the critical values for h can be calculated 
for each tube. 

The experimental results are displayed in Figure 1.8. There is 
excellent agreement except at large values of B, where the theory under 
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under-predicts the critical stress. One expects this, because at these 
large values of B the puddle itself nearly fills the tube before slugging 
occurs. As seen in Figure 1.4(d), when the puddle level is near the top, 
the transverse radius of curvature R£ is on the average much smaller than 
the tube radius Rr that is used to approximate it in the theory. 

The slugging criterion for "Dee-shaped" vapor spaces was measured 
experimentally for shapes corresponding to a typical ratio of the wick 
thickness to tube diameter of 0.465. Long stainless-steel strips of the 
appropriate width were inserted into each glass tube and held in place 
with a spring. Slugging of the resulting "Dee-shaped" vapor space was 
measured with the "Dee" on the top and the bottom, which corresponds to 
a horizontal wick, and on the side, which corresponds to a vertical wick. 
The results are displayed in Figure 1.10. In the case of the "Dee," the 
Bond number is based on one-half of the hydraulic diameter D|^ instead of 
the tube radius. The faired lines through the data of Figure 1.9 are 
given by an empirical data fit 

,t/d, = ^ ^ (1 - 

where the empirical constants are given in Table 1.1. These empirical 

(4) 

results are incorporated in the computer program GRADE' ' so that optimum 
graded-porosity wicks can be designed that just prevent the occurrence of 
a liquid slug. 

Table 1.1. Constants for the Empirical Curve Fit for 
the Slugging Condition for "Dee-Shaped" 

Vapor Spaces 

hj/Dh = 1/B + (a/B^)(l - n 




i 

X 

11 

Dee on top 

.820 

.40 

1.5 

2.9 

Dee on bottom 

.68 

.40 

1.5 

1.9 

Dee on side 

1.5 

.5 

1.0 

1.2 
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a / (a + H) 
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9 


Figure 1.8. Slugging Criterion for Closed-Ended Cylindrical Tubes 



26263-6025-RU-00 



Figure 1.9. Schematic Diagram of Apparatus to Measure 
the Slugging Criterion for Closed-Ended 
Cylinders. Slug is shown formed in the 
smallest tube. 
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Figure 1.10. Empirical Sluaging Criterion for "Dee-Shaped" Vapor Spaces 
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2.0 THE CAPACITY OF CIRCUMFERENTIAL GROOVES 

Many heat pipes use a wick structure for the axial transport of 
liquid, and circumferential grooves for the transport across the inner 
tube wall in the evaporator and condenser regions. This combination 
provides both high axial transport and heat-transfer coefficients. 

Our particular interest is the use of circumferential grooves with the 
graded-porosity slab wick; however, the results of this section apply 
to any heat pipe that uses circumferential grooves. 

A sound mathematical model of circumferential grooves is required 
(4) 

for the program GRADE' ' . which designs graded-porosity wicks, because 
the evaporator grooves are the limiting factor for the transport capacity. 
GRADE repeatedly integrates the governing differential equations. Each 
integration uses an assumed value for the axial heat- transport rate, and 
a subroutine DRY reports whether the circumferential grooves fail. If 
they do, the assumed heat-transport rate is reduced, and if they don't 
the rate is increased. Then the integration is repeated. A binary 
search on the rate is used to design a graded-porosity wick with the 
highest axial heat- transport capacity that will not fail the circum- 
ferential grooves. 

In this section, we present the mathematical model for the sub- 
routine DRY. The method used is similar to that used by Gier and 
f5l 

Edwards' s except that here the more general trapezoidal groove rather 
than the triangular groove is considered, and a general wick structure 
is considered rather than one that feeds the grooves only at the bottom. 

Figure 2.1 depicts the geometry under consideration. To be 
specific, a vertical slab wick is shown. However, the method applies 
to any wick configuration. We focus attention on a single region of 
grooves that is fed liquid from the wick at angular locations ej and 
02 that are measured from the top (positive in the counterclockwise 
direction). The total heat input per unit length along the heat pipe 
is Q, and a fraction Fq is assumed to enter between ej and 02 . 
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Figure 2.1. Liquid Enters a Groove s . the Rate M, 
at Location ar ’ it tF Rate M 2 at 
Location 02 . It is uniformly evaporated 
between ei and C2. 
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We assume that the heat Is evaporated uniformly along the groove, 
thus the amount of liquid evaporated between location e and e+de 1 s 

dm(e) = -(Fp QAN) [de/( 02 - 6 i)] (2.1) 

where m(9) Is the local mass flow through the groove assumed positive 
In the counterclockwise direction, N Is the number of grooves per unit 
length and X is the latent heat. Integration of Eq. (2.1) over the 
Interval (e^ , 62 ) yields 

(mj + m 2 ) = Fp Q/XN (2.2) 

where mj and m 2 are the absolute values of the flow rates Into the 
groove at and 62 - One location between Oj and 62 Is a stagnation 
point where the flow rate Is zero. By integrating Eq. (2.1) from &i and 
then 62 to 6 g, we obtain the result 

mi/m2 = (6^-di}/^e2-e^) (2.3) 

Unless the region of grooves under consideration is situated symmetrically 
in a gravitational field or they are operating in zero gravity, in which 
case mj equals m 2 , then one does not know a priori the values of m^ and 
m 2 . Their sum is given by Eq. (2.2) and the unknown location of the 
stagnation point fixes their ratio by Eq. (2.3). Combining these equations 
we obtain expressions for m^ and m 2 in terms of e^: 


= (Fp Q/XN) (05-0i)/(02-0i), 

(2.4) 

“ (Fp Q/XN) (02“05)/(02”®i) • 

(2.5) 


Thus our lack of knowledge of the values of mj and m 2 is shifted to the 
unknown location of the stagnation point. Integrating Eq. (2.1) from 
6 , with the use of Eq. (2.4) we obtain for the mass flow rate 

m(e) = (Fp Q/xN) (e5-0)/(e2-0i) (2.6) 
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As liquid flows through the groove, its pressure drops relative to 
that of the vapor. The pressure difference is sustained by the curvature 
of the free-surface meniscus. The rate of increase of the liquid-vapor 
pressure difference P is governed by the expression 


^ = - Ry ^ (1/2 p u 2)/D^ - pgR^ sin e 
e 


(2.7) 


where 


In terms 


Rj is the inside radius of the tube 
K is the friction-factor coefficient 
Rg is the Reynolds number based on the hydraulic diameter 
p is the liquid density 

U is the average liquid velocity, 

of the mass flow rate given by (2.6), we write 


dP _ 

?e ■ 


Rj K V 


(Fq Q/xN)(eg-e)/( 02 - 0 i) - pg Rj sin 0 


( 2 . 8 ) 


This equation must be integrated numerically because K, Dj^, and A vary 
depending on the shape of the meniscus, which in turn depends on P. 

Since typically the Bond number is very small, we assume that the meniscus 
is circular, and thus the liquid in the grooves is as shown in Fig. 2.2. 
The radius is given in terms of the local l!quid-vapor pressure differ- 
ence P and surface tension o by: 


R ' o/P (2.9) 

When the vapor-liquid pressure difference is low, the meniscus is attached 
to the top of the groove. As the pressure difference increases, the radius 
of curvature of the meniscus decreases. When the contact angle x reaches 
the wetting angle o, a further increase in the pressure difference results 
in the meniscus receding into the groove as shown in Fig. 2.2(b). Groove 
failure is taken to occur when the meniscus reaches the groove bottom. 
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Figure 2.2. Cross Section of a Groove 

(a) Meniscus of liquid attached to top 

of groove with contact angle x. I 

(b) Meniscus receded into groove {contact 
X equals wetting angle a). 
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The area A and hydraulic diameter are easily calculated for 
a given value of P once R is calculated from Eq. (2.9). For the 
calculation of the friction factor coefficient, we rely on the work by 
Hsu^^^ who used numerical techniques to calculate the friction factor 
for a circular meniscus in a trapezoidal groove. The subroutine DRY 
contains a subroutine FACTOR that is devoted to producing an empirical 
curve fit to Hsu's numerical results. 

The boundary conditions for Eq. (2.8) are the values of the vapor- 
liquid pressure difference at the two feed locations, which are related 
hydrostatically to the stress S at the tube center; 

Pi = S + pgRj cos 01 , 

P 2 = S + pgRj cos 02 , (2.10) 

Equation (2.8) with the boundary conditions (2.10) constitute a two- 
point boundary-value problem where the unknown stagnation point 0 ^ is an 
eigenvalue. One can attempt a solution by repeatedly integrating from 
01 to 02 . each time adjusting the value of 0 ^ until both boundary condi- 
tions are met. A solution will be possible if the heat input Q is below 
the critical value at which the groove fails. Such an iterative search 
for a solution is time consuming, especially when the computer subroutine 
that would carry it out is itself called repeatedly. It is often possible 
to decide whether or not Q is above the critical value based on a single 
integration. Consider a curve constructed as follows. We assume a 
value for 0 ^ and then integrate F.q. (2.8) from 0i to 0^. In Fig. 2.3 we 
plot the vapor-liquid pressure difference P^ at the stagnation point as 
a function of 0 . The curve ends when P is sufficiently high to fail 
the groove. Similarly, we integrate from 02 to 0^ and again plot P^ as 
a function of 0 ^. If Q is below the critical value for groove failure, 
the two curves will intersect at the actual stagnation point, as shown 
in Fig. 2.3(a). If, on the other hand, Q is above the critical value, 
they will not intersect as shown in Fig. 2.3(b). 

We now can describe the procedure used to assess whether or not the 
groove has failed. A trial value of is assumed half way between 0 i and 
02 , and we attempt to integrate to it first from and then from 62 . If 
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U4 



01 02 


ASSUMED STAGNATION POINT LOCATION 0^ 


UJ 



©1 ©2 


ASSUMED STAGNATION POINT LOCATION © 


Figure 2 . 3 . Plots of the Vapor-Liquid Pressure Differ- 
ence at an Assumed Stagnation Poir.t Location 
©5 by Integrating from Both 0-| and ©2 to ©5. 
(a) Case where grooves function, (L) Case 
where they dry out. 
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two values of are successful lly calculated without groove failure, then 
evon if they are not tnc ' me ve havi^ ated w..-it the two cut >/es 

are as in Fig. 2.3(a) and they must intersect. Therefore, we conclude that 
the groove has not failed. If both attempts to integrate to fail, then 
we have demonstrated that the two curves are as in Fig. 2.3(b) and they 
do not intersect. Therefore, we conclude th-.c the groove has failed. 

Only when integration to the trial value for is successful from one 
direction but not from the other, can we draw no conclusion. When this 
happens, a new trial value for is selected, and the procedure is repeated 
until a conclusion can be drawn. Usually only a few trys are needed at 
most. 

The procedure just described is incorporated in the computer sub- 
routine DRY, which reports whether the grooves have failed for a given 
heat input and capillary stress at the tube center. Sample results are 
shown in Fig. 2.4. The primary input to DRY, besides parameters describ- 
ing the groove geometry and fluid properties, is the heat input per unit 
length of heat pipe and the capillary stress, which is specified by 
the radius of curvature R^^p of the free surface in a groove under the 
capillary stress S: 


S = 


o/R 


cap' 


The region above the curve in Fig. 2.4 corresponds t .lues of R^^p and 
Q|^ where the grooves are failed, while the region btlow it corresponds to 
values where the grooves are not failed. 
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(10’^ M) 


Figure 2.4. Sample of Groove Dryout Calculations 


- 30 - 



26263-6025-RU-00 


3.0 VAPOR FLOW IN HEAT PIPES 

We consider in this section the effect of vapor flow on wick- 
limited heat pipes. Vapor flow is often an important design considera- 
tion since the vapor-liquid pressure difference that must be supported 
by the capillary pressure of the work is the sum of the pressure drop 
in the vapor as well as in the liquid. If a heat pipe has a negligible 
vapor-flow pressure drop at its maximum heat-transport rate, then it is 
often possible to increase its capacity by increasing the cross-sectional 
area of the wick. In the limit of this process where the wick completely 
fills the heat-pipe cross section, the capacity, of course, falls to zero. 
Therefore, an optimum cross-sectional area for the wick exists and when 
other factors such as radial heat transfer and vapor-space slugging are 
not a consideration, this optimum occurs when the vapor-flow and liquid- 
flow p.essure drops are of the same magnitude. An accurate calculation 
of the vapor-pressure profile is especially important for a graded- 
porosity wick because it is designed to just sustain the vapor-liquid 
pressure difference along its length rather than only at the evaporator 
end. 

The vapor-flow calculation is difficult because heat pipes often 
operate in a regime where there is an essential interaction between 
viscous and inertial forces. This has been the subject of considerable 
study which has been reported in the literature. Bankston and Smith^^\ 

/o\ 

and Tien and Rohani' for example, resorted to a numerical solution 
of the full Navier-Stokes equations, while Busse^^^, Galowin et 
Quail and Levy ^^^,13) Gupta and Levy^^^^ based their analysis on 
the boundary-layer equations. Their solutions are limited to either 
circular, plane-parallel or annular vapor-space cross sections. We are 
primarily interested in, rather, the "Dee-shaped" vapor space formed by 
a graded-porosity slab wick in a circular tube. Also, the mathematical 
complexity of the solution of either the Navier-Stokes or the boundary- 
layer equations makes such an approach unsuitable for routine computer- 
design calculations. In section 3.1 we present a simplified approximate 
mathematical model that is based on the results of earlier more exact 
theories. It has the virtue of providing a fast calculation of the 
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static-pressure profile of the vapor in a Dee-shaped vapor space while 
retaining the effect of the inertial -viscous interaction. In section 
3.2, we present the results of an experiment that actually measures the 
pressure profile in the vapor. 

3.1 THEORETICAL MODEL FOR VAPOR FLOW IN A DEE-SHAPED SPACE 


3.1.1 Flow in the Adiabatic Section 


A rigorous analysis of the actual Dee shape would require at least 
a three-dimensional solution of boundary- layer equations, which is far 
beyond the scope of a design-oriented subroutine. Therefore, as our 
first step we approximate the Dee shape by a two-dimensional configuration 
that consists of two flat parallel plates. One is always adiabatic and 
models the surface of the wick. The other, which models the tube wall, 
is either a surface of evaporation, condensation, or it is also adiabatic. 
A question arises as to the plate separation that will best approximate 
a Dee-shape. The criterion that is used requires that the flow in both 
the actual Dee-shaped passage and between the parallel walls have the 
same average velocity and the same pressure gradient dP/dx in adiabatic 
regions of fully developed flow. If we let the subscript ( )p refer to 
the Dee-shape and ( )j| refer to the parallel plates, then the criterion 
for the plate spacing is 


dP 

dx 


_ 96 
■ 1 ??. 


1/2 p u2 

,~U77" 


^^D* *^®D^ 1/2 p u^ 




(3.1) 


where fp is the friction factor. Re is the Reynolds number based on the 
hydraulic diameter D, p is the vapor density and u is the average 
velocity. After we have substituted twice the plate spacing b for the 
hydraulic diameter of parallel plates, equation (3.1) reduces to 

b = 2 /6/fjj Rejj Dp (3.2) 
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The friction factor fjj depends on the ratio of wick thickness to 
tube diameter, and In general It must be found from the numerical 
solution of Poisson's equation: 


3^U 4 . 3^U 


(3.3) 


where here the z axis is taken as the flow direction, x and y are trans- 
verse coordinates, u Is the flow velocity, p Is the pressure and y Is 
the viscosity. Equation (3.3) Is written In dimensionless form by 
Introducing the dimensionless variables 


U 


z y 


(3.4) 


x/d. 

(3.5) 

y/d. 

(3.6) 


where d 1s the tube diameter, and we obtain. 


^ = .1. (3.7) 

3X2 3Y2 

The boundary conditions require U to vanish on the boundary. Once we 
obtain a solution, the friction factor is related to the average 
dimensionless velocity by 

¥«D " . (3.8) 


Equation 3.7 is written In a finite-difference form by overlaying the 
flow region on an N x N rectangular grid. By symmetry, only half of 
the region must be used. Solution for the velocity U,.. at the (1,j) 

* J 

node Is obtained by Gauss-Seidel Iteration of the finite difference 
approximation to Eq. (3.7). The resulting friction-factor Is shown 
in Figure 3.1. 

Analytical results are possible for w/d = 0 and asymptotically for 
w/d-^-l. The result of Oon = 62.8 for w/d = 0 Is a special case of an 
analysis by Sparrow' An asympt^'c solution Is possible in whe 
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Figure 3.1 Friction-Factor For A "Dee-Shaped" Conduit 
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limit w/d-^) because In Eq. (3.7) can be neglected. The result- 

ing friction factor 

i/M V'O = 'W’ ' “-2 <3-«) 

The numerical solution in Fig. 3.1 does not quite agree with the analytical 
results for w/d = 0 and 1. This is attributed to the approximation of 
the exact .differential equation by the finite difference equation. In 
fact, the friction factors shown in Fig. 3.1, which were calculated for 
a 30 X 30 mesh, increased somewhat as a finer mesh was used. 

We were surprised to find that the friction-factor is practically 
independent of w/d. For practical calculations, the average of the 
analytical results for w/d = 0 and w/d = 1 is recommended. This value is 
62.5/Rejj. 

3.1.2 Flow in Condensing Sections 

In the condensing section the continual loss of mass flow to the 
condenser wall requires a deceleration of the main flow which is brought 
about by a positive pressure gradient in the direction of flow. The 
positive pressure gradient tends to cause a flow reversal near the wall, 
where viscosl^ iS'ings the tangential velocity to a low value so that 
inertia effects are overcome by the positive pressure gradient. This 
behavior is particularly true near a noncondensing surface, such as that 
of the slab wick. The tendency toward flow reversal lowers or indeed 
reverses the viscous shear stress which otherwise acts against the flow. 

Consider an element dz long of a heat pipe with walls forming a 
cylinder in the general mathematical sense, i.e., parallel wall elements 
but not necessarily with a circular cross-section. A force-momentum 
balance requires the following: 

^ p <1 Ac = - ^ jT P d Ac - (3.10) 

where A is cross-sectional area, is wall shear stress and C is 
perimeter. The area- integrated terms will be written as overscored 
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quantities as follows 


l 

l 


P dA 
p dA = 



pA 


A 


(3.n) 

(3.12) 


Additionally, e. velocity-squared shaped factor Is Introduced 


Fg = ^/i"2 


(3.13) 


A circumferentially averaged friction factor f (not to be confused with 
friction-factor coefficient to follow) Is defined: 


/’w 


dC - C •g- pu^ 


(3.14) 


Hydraulic diameter Is also Introduced 

D = 4 A/C (3.15) 

With these definitions Eq. (3.10) takes the form 

Equation (3.16) can be Integrated over an entire length of a condenser 
with closed end at x = L to yield 

Pl “ Po “ %,o '"“o ■ j[ ^ 1 

For a condenser with a uniform rate of condensation 

u(z) = Uq [1 - 2 /L] (3.18) 

Hence, Eq. (3.17) becomes 

Pl ■ Po ■ %.o P®' - W P"o I * 
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where 


i - x/L. 


A pressure-recovery coefficient is defined by 


r ' (Pl-P,)/(j ou^) 


Re . a . 

V 


(3.20) 


(3.21) 


(3.22) 


and a friction-factor coefficient K is defined by 

f = K/Re 
where Re is Reynolds number 

in 

(3.23) 

ai^d V is the kinematic viscosity. Equation (3.19) can thus be written 

d5 (3.24) 


^ %,o ‘ D R^ j[ 


Finally an average value of friction factor coefficient K can be defined by 


K = 2 


/’ K[1-0 


dC 


(3.25) 


Equation (3.24) can then be written 


r = 2 F 


L X 

s ,0 Re^ 


(3.26) 


Inlet Reynolds number Re^ and radial Reynolds number 


Re. 


u 0 
w 


r V 

may be used interchangeably, because 


% \ CL = 4 u^ L/D 


Re^ = 4(L/D) Re^ 


(3.27) 

(3.28) 

(3.29) 
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Hence, Eq. (3.26) can be written 

^ %,o ‘ rli; 

or 

R • 8 Re^ [2 ,, - r] 

It is clear that, given a value of pressure recovery coefficient 
r from the literature, as a function of Re^ or Re^ alone for a complete 
condenser, and given the inlet velocity distribution so that Eq. (3.13) 
can be used to fix ^ Eq. (3.31) can be used to find the average 
friction factor coefficient K . It is suggested that K be used in 
approximate local pressure-distribution calculations merely by employing 
it in Eq. (3.16) in place of the true K, as follows: 

i = ■ %,o ^ - fe ^ (3-32) 

A constant value of is also used for simplicity, which is taken as 
the value for fully developed channel flow. While only an approxi- 
mation, Eq. (3.32) will integrate to yield the correct result for Pj^, 
since that result is built in. It is recommended that calculations be 
made with the approximate Eq. (3.32) and compared to more exact numerical 
or experimental results. 

Gupta and Levy^^^^ report experimental results obtained with 
silicone oil flowing between a solid outer cylinder and a porous inner 
one. Results were obtained for an inner to outer radius ratio of 0.83. 
Due, presumably, to flow separation and turbulence, r achieved a value 
of only + 0.5 as radial Reynolds number Re^ becomes large. The form of 
Eq. (3.30) and the behavior displayed in Figures 10 and 11 of Ref. (14), 
suggested fitting the results with the following form of equation 

r = A - ^ exp {-[a Re^ + b Re/]} (3.33) 

r 


(3.30) 

(3.31) 
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In the limit of large Re^ this relation goes to the constant A (equal 

to 0.5 from the experimental results). At small values of Re^, r goes 

to -B/Re^ so that Eq. (3.31) gives K = 8B; hence B must be the ordinary 

Poiseuille value of K divided by 8; B=12 for parallel plates and 8 for 

the circular tube. The constants a and b govern the rapidity with 

which the relation goes from the -B/Re„ asymptote to the constant A 

(14) ^ 

asymptote. Gupta and Levy' ' also report calculated values for 
parallel plates. 

Table 3.1 shows their values and the fit according lo Eq. (3.33). 

The fit is within the error in reading the curve in Ref. (14). Table 3.1 
also shows the values of K corresponding to the values of r according to 
Eq. (3.31). The values show that, due to the fact 2F is considerably 

o 9 (j 

larger than A, K increases strongly with increasing Re^. 

3.1.3 Flow in the Evaporating Section 

In the evaporating section of a heat pipe the continual addition of 
mass flow requires an acceleration of the main flow which is brought 
about by a negative pressure gradient in the direction of flow. Viscous 
drag by the wall increases the magnitude of the pressure gradient. At 
very high rates of evaporation this factor becomes unimportant, as the 
velocity gradient at the wall tends to zero. At low rates of evaporation 
the viscous drag governs. 

The analysis proceeds in the same manner as in the case of condensing 
regions to the approximate governing equation 

t = - fs.L i - Ir # (3.34) 


Where Fj ^ 
given by 


is the value of F^ at the exit of the evaporator and K is 



(3.35) 


In this case, the subscript ( )o refers to the closed end of the 
evaporator and ( )|^ refers to the exit. The overall pressure drop 
coefficient defined by 
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TABLE 3,1 

Pressure Recovery Coefficient 
Versus Radial Reynolds Number 
For Flat Plates 


Radial 
Reynol ds 

Recovery Coefficient r 

Average Friction 
Factor Coefficient 

Number 

Re^ 

■ 

1 

Eq. (3.33)** 

K (Eq. 3.31)* 

0 

N.A. 

-12/Re^ 

96 

1 

N.A. 

-9.9 

98 

2 

N.A. 

-3.9 

100 

3 

N.A. 

-1.92 

104 

4 

-0.9 

-0.97 

108 

5 

-0.4 

-0.43 

113 

6 

-0.1 

-0.09 

119 

7 

0.1 

0.13 

127 

8 

0.3 

0.28 

136 

10 

0.5 

0.45 

156 

12 

N.A. 

0.53 

179 

14 

N.A. 

0.57 

205 

16 

N.A. 

0.59 

232 

00 


0.60 

14.4 Re^ 



6/5 


** A = 0.6, B - 12, a = 0.13, b = 0.008 


- 40 - 










26263-6025-RU-00 


r = (Pq-Pl)/?- P\ 
is related to K and ^ by 

K = 8 Re^ [r - 2 Fj ^3 

We again fit Gupta and Levy*s^^^^ results this time with 

r = A + exp {- [a Re^ + 6 Re^ 
r 

As in the case of the condensing region, B is taken as the Pouiseuille 
value of the friction-factor coefficient (96 for flat plates) divided 
by 8 in order to recover Pouiseuille flow in the limit of small Re„. 
Unlike the condensing region, however, in the evaporator at large values 
of Re^ the viscous term should vanish, thus we take 

A = 2 Fj L (3.39) 

Table 3.2 shows values of r from Ref. (14), the fit of Eq, (3.38), and 
K according to Eq. (3.37). 

3.2 EXPERIMENTAL VERIFICATION OF THE VAPOR-FLOW MODEL 

In order to verify the theoretical model of vapor flow that is 
incorporated in GRADE II, a task was undertaken to experimentally 
measure the vapor pressure profile in a heat pipe that closely resembles 
ones typical of a spacecraft application. Figure 3.2 depicts the cross 
section of the experimental heat pipe. It is 188cm (74 in.) long with 
19 pressure taps spaced 10.16 cm (4 in.) apart starting one inch from 
the evaporator end. 

The pressure taps present a unique design challenge. The usual 
practice in measuring the static pressure of a gas stream at a wall is 
to drill a small -diameter hole on the order of 0.C25 cm (0.010 in.) in 
the wall, which is connected to a pressure transducer by means of tubing. 
Such a small hole in a heat-pipe wall, however, would fill with liquid 
and the resulting capillary pressure would influence the measurement. 


(3.37) 

the equation 
]} (3.38) 
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TABLE 3.2 


Pressure Drop Coefficient and 
Average Friction Factor Coefficient 
Versus Radial Reynolds Number 


Radial 

Reynolds 

Number 

Pressure Drop Coefficient 

— 

Average 

Friction Factor 
Coefficient** 

Gupta & Levy [14] 

Eq. (3.38)* 

0 


12/Re^ 

96 

2 

8 

8.21 

78.6 

3 

6.4 

6.26 

71.1 

4 

5.5 

5.31 

64.4 

5 

4.9 

4.76 

58.2 

6 

4.4 

4.40 

52.7 

7 

4.1 

4.15 

47.7 

8 

3.9 

3.97 

43.1 

10 

3.7 

3.74 

35.3 

12 

3.6 

3.60 

28.9 

14 

3.5 

3.51 

23.7 

16 

3.4 

3.45 

19.4 


* A = 3.3, B =• 12, a = 0.10, b = 0 

** Based upon F5 l ~ 
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Figure 3.2 


Cross Section Through Pn 
Pipe for Measuring Vapor 
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As shown in Fig. 3.2, our solution to the problem is to use a large- 
diameter pressure tap whose capillary pressure is insufficient to support 
a liquid slug. The mathematical criterion in terms of the surface 
tension a, liquid density p, pressure tap I.D., d, and heat pipe I.D., 

D, is 

4o/d < pgD. 

(We have assumed that the heat pipe is not overfilled, that is, there 
is no puddle along the bottom.) The wall hole is covered on the inside 
by a patch of 0.00127 cm (0.0005 in.) foil with an 0.010 inch diameter 
hole in it. The hole stays clear of liquid by menisci coalescence, 
which is the same mechanism on which the priming foil is based. A 
sample pressure tap was fabricated and tested. A liquid slug inten- 
tionally introduced into the tube was observed to drain away and the 
foil hole cleared itself by menisci coalescence. 

The pressures at each tap were measured relative to the first 
pressure tap at the evaporator end with a Validyne variable-reluctance 
differential pressure transducer (0.1 psid full scale). One side of 
the transducer was connected to the evaporator-end tap, and the other 
side was connected to each of the other taps via a valve. A differ- 
ential pressure measurement is made to a particular pressure tap with 
only the va^ve to that tap open. 

The experimental results and the corresponding theoretical cal- 
culations by GRADE II are displayed in Fig. 3.3. MethanoT is used as 
the working fluid since its low vapor density at room temrerature 
results in a high vapor velocity and hence a relatively high pressure 
drop. Figure 3.3 shows good agreement between theory and experiment 
at a 20 watt heat-transfer rate and 10®C vapor temperature. The radial 
Reynolds number of 0.256 in the condenser section and 0.235 in the 
evaporator section shows that the flow is dominated by viscous effects. 

To increase the inertial effect, the heat transfer rate was increased 
and the vapor temperature was decreased. The results for 30 watts and 
10“C shown in Fig. 3.3 was the highest radial Reynolds number that was 
obtained before the pressure readings in the condenser section became 
unstable and large fluctuations occurred. 
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Instability of the flow field in the condenser region has been 
measured previously^^^’^^\ but at higher radial Reynolds numbers. 

We cannot completely account for the instability at such low Reynolds 
numbers. One factor that could strongly contribute, however, is the 
fact that the vapor-space surface is very rough, especially the slab 
wick surface. Such roughness would be expected to induce unsteady 
flow at a much lower Reynolds number. 
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